The properties of D and D* mesons in the nuclear medium 
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Abstract 

We study the properties of D and D* mesons in nuclear matter within a simultaneous self-consistent 
coupled-channel unitary approach that implements heavy-quark symmetry. The in-medium solution ac- 
counts for Pauli blocking effects, and for the D and D* self-energies in a self-consistent manner. We pay a 
special attention to the renormalization of the intermediate propagators in the medium beyond the usual 
cutoff scheme. We analyze the behavior in the nuclear medium of the rich spectrum of dynamically-generated 
baryonic resonances in the C = 1 and S = sector, and their influence in the self-energy and, hence, the 
spectral function of the D and D* mesons. The D meson quasiparticle peak mixes with E C (2823)7V~ 1 and 
S C (2868)7V _1 states while the Ac(2595)iV _1 mode is present in the low -energy tail of the spectral function. 
The D* spectral function incorporates J = 3/2 resonances, and £ C (2902)7V~ 1 and A C (2941)7V -1 fully de- 
termine the behavior of D* meson spectral function at the quasiparticle peak. As density increases, these 
resonant-hole modes tend to smear out and the spectral functions get broad. We also obtain the D and 
D* scattering lengths, and optical potentials for different density regimes. The D meson potential stays 
attractive while the D* meson one is repulsive with increasing densities up to twice that of the normal 
nuclear matter. Compared to previous in-medium SU(4) models, we obtain similar values for the real part 
of the D meson potential but much smaller imaginary parts. This result can have important implications 
for the observation of D°-nucleus bound states. 



PACS numbers: 11. 10. St, 12.38. Lg, 14.20.Lq, 14.40.Lb, 21.65.-f 



I. INTRODUCTION 



The interest on the properties of open and hidden charmed mesons was triggered more than 20 
years ago in the context of relativistic nucleus-nucleus collisions in connection to the charmonium 
suppression [1] as a probe for the formation of Quark-Gluon Plasma (QGP). The experimental 
programme in hadronic physics of the future FAIR facility at GSI 0] will move from the light 
quark sector to the heavy one and will face new challenges where charm plays a dominant role. 
In particular, a large part of the PANDA physics programme will be devoted to charmonium 
spectroscopy. Moreover, the CBM experiment will extend the GSI programme for in-medium 
modification of hadrons in the light quark sector, and provide first insight into the charm- nucleus 
interaction. 

The primary theoretical effort is to understand the interaction between hadrons with the charm 
degree of freedom. Charmed baryonic resonances have received recently a lot of attention motivated 
by the discovery of quite a few new states by the CLEO, Belle and BABAR collaborations [3, 4, 
SUSS]- Whether those resonances have the usual qqq structure or qualify better as being 
dynamically generated via meson-baryon scattering processes is a matter of strong interest. In 
fact, the unitarization, in coupled-channels, of the chiral perturbation amplitudes for scattering of 
CP octet Goldstone bosons off baryons of the nucleon l/2 + octet has proven to be quite successful 
in the charmless sector @, EI EI HI EI EI EI EI El, EI EI M, Q 0, 0, 0, 0, 0, 0]- 
The modification of the various meson-baryon amplitudes for the case of finite temperature and / or 



nuclear density has also attracted a lot of attention and has been carefully discussed [23, |30j, 131 



3J, |33j,|34|. 



The extension to the charm sector of the unitarized meson-baryon method was attempted in a 



first exploratory work in Ref. 35], where the free space amplitudes were constructed from a set of 



separable coupled-channel interactions obtained from chirally motivated lagrangians upon replacing 
the s quark by the c quark. A different approach resulting from the scattering of Goldstone bosons 
off the ground state l/2 + charmed baryons was pursued in , but the substantial improvement in 
constructing the meson-baryon interaction in the charm sector came from exploiting the universal 
vector- meson coupling hypothesis to break the SU(4) symmetry [36]. The t-channel exchange of 
vector mesons (TVME) between pseudoscalar mesons and baryons preserved chiral symmetry in 
the light meson sector keeping the Weinberg- Tomozawa (WT) type of interaction. An extension 
to (f-wave J = 3/2~ resonances was developed in [37l ] . while some modifications over the model of 
Ref. [3^ | were implemented in Ref. 38], both in the kernel and in the renormalization scheme. More 
recently, there have been attempts to construct the DN and DN interaction by incor por ating the 
charm degree of freedom in the SU(3) meson-exchange model of the Jiilich group 

Nuclear medium modifications were then incorporated in order to study the properties of 
charmed mesons in nuclear matter, and the influence of those modifications in the charmonium 
production rhythm at finite baryon densities. Possible variation of this rhythm might indicate 
the formation of the QGP phase of QCD at high densities. Previous works based on mean-field 
approaches provided important mass shifts for the D and D meson masses 41 , 0, [4I 44 ] . Some 
of those models have been recently revised 45, 46j ]. However, the spectral features of D mesons in 
symmetric nuclear matter were obtained for the first time in the exploratory work of Ref. 35J] , while 
finite temperature effects were incorporated later on in Ref. [43]. Afterwards, within the SU(4) 
TVME model of Ref. 36], the properties of the {D,D) and (D s , D s ) mesons were analyzed in 1481 1. 
and in 
Ref. 



38l . |43]. In this latter reference, the kernel and the renormalization scheme employed in 



361 ] were modified. 



However, those SU(4) TVME inspired models are not consistent with heavy-quark symmetry 
(HQS), which is a proper QCD spin-flavor symmetry that appears when the quark masses, such 
as the charm mass, become larger than the typical confinement scale. As a consequence of this 



2 



symmetry, the spin interactions vanish for infinitely massive quarks. Thus, heavy hadrons come 
in doublets (if the spin of the light degrees of freedom is not zero) , which are degenerated in the 
infinite quark-mass limit. And this is the case for the D meson and its vector partner, the D* 
meson. 

In fact, the incorporation of vector mesons into the coupled-channel picture has been pursued 
very recently in the strange sector. On one hand, vector mesons have been incorporated within 
the hidden-gauge formalism. Within this scheme, a broad spectrum of new resonant meson-baryon 
states have been generated 5(], 51, Hsf. O n the other hand, the WT meson-baryon chiral Lagrangian 
has been also extended to account for vector meson degrees of freedom by means of a scheme 
that starts from a SU(6) spin-flavor symmetry Lagrangian and that incorporates som e sy mmetry 



breaking corrections determined by physical masses and meson decay constants [53J, [54|, [55|, |56( |. 
The corresponding Bethe-Salpeter equation reproduces the previous SU(3)-flavor WT results for 
the lowest-lying s- and d-wave, negative parity baryon resonances and gives new information on 
more massive states, as for example the A(1800) or A(2325) resonances. The extension of this 
scheme to four flavors, incorporating the charm degree of freedom, was carried out in Ref. (5?J 
and it automatically incorporates HQS in the charm sector improving in this respect on the SU(4) 
TVME models, since D and D* mesons are thus consistently treated. One of the distinctive 
differences of this approach with respect to those that built in the SU(4) TVME model can be 
drawn in the wave function content of the resonances. Thus, for instance, the dynamics of the 
lowest lying resonance A c (2595) is completely dominated by the DN channel in the SU(4) TVME 
model of Ref. (3^1, while it turns out be largely a D*N state within the SU(8) scheme of Ref. 57]. 
Such differences might have also a direct influence on the dynamics at finite densities. 

Thus in this work, we aim to investigate the nuclear medium effects in hadronic systems with 
charm one (C = 1) and no strangeness (S = 0) within the SU(8) model derived in [57j. In 
particular, we study the dynamically-generated baryonic resonances in the free space as well as in 
the nuclear medium in order to analyze how the masses and widths are modified with density. We 
also study the D and D* mesons self-energies in the nuclear medium, calculating their spectral 
functions for a variety of densities and the corresponding optical potentials. A novelty of our work 
is that we simultaneously obtain, in a self-consistent manner, the D and D* meson self-energies. We 
then compare our results with the previous ones obtained within SU(4) TVME schemes and other 



simple models 



491 ]. paying also an special attention to the regularization of the 



more 

intermediate propagators in the medium beyond the cutoff method. The nuclear medium effects for 
hadronic scattering amplitudes and for hadrons propagators are of interest for the understanding 
and correct interpretation of the data obtained in heavy-ion collisions where the high nuclear 
densities reached can substantially change the properties of the involved hadrons. 

To end this introduction, and before we start deriving the D and D* self-energies inside of a 
nuclear medium, we would like to make a general reflexion to better situate this work. Though the 
model of Ref. 57J is possibly the best one existing in the literature for describing the free space 
C = 1, S = meson-baryon elastic scattering at low energies 1 , it does not provide, as it is the 
case for the rest of the available models, the correct analytical properties of the scattering matrix, 
including s, t and u— cuts and proper crossing symmetry. Only the real theory, presumably QCD, 
could do that. However, one should bear in mind that deficiencies of the free space model would 
definitely affect to results that will be presented here for amplitudes embedded in cold nuclear 
matter. Since we aim at describing resonances, it is important to use a model consistent with the 
unitarity cut as that of Ref. [57J], while for rest of analytical properties one hopes to be partially 



As outlined above, it provides a scheme for four flavors and for pseudoscalar and vector mesons which reduces to 
the WT interaction in the sector where Goldstone bosons are involved, and that incorporates HQS in the charm 
sector. 
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taken into account thanks to low energy constants (subtraction constants...). To ignore a cut that 
would cover totally /partially the studied resonance region is not so much important 2 , and what 
is more relevant is the proximity of a resonance to a branching point, where the amplitude might 
vary more rapidly, together with the strength of the coupling of the resonance to the given channel. 
The free space model of Ref. [53] only properly accounts for the two body unitarity cut, and it 
thus suffers from some limitations. As a matter of example, this model does not account for the 
dynamics of an intermediate nucleon in the u— channel diagram DN — > A c p, involving two p— wave 
couplings (DNA C and NNp), that gives rise to a cut extending from 1.55 to 2.74 GeV. The lowest 
energy branch point should not suppose a serious problem. On the other hand, A c p selects isospin 
1=1, and if one looks at Tables VI and IX of Ref. [57], one sees no resonances excessively close to 
the branch point 2.74 GeV, with a dominant coupling to that channel. Thus, though some effects 
will exist, we do not expect drastic changes that could lead us to think that the interaction model 
of Ref. [53] is unrealistic, specially close to threshold where the p— wave couplings will be negligible. 

Other contributions no considered in the model of Ref. [57J are those driven by the exchange 
of a pion in the t— channel. Among all of them that contributing to the DN —* D*N amplitude 
could be quite relevant, since it might induce a pronounced energy dependence, difficult to account 
by means of the low energy constants of the model 3 . This contribution involves again two p— wave 
vertices, which would vanish at the D*N threshold. Below threshold, its contribution would 
heavily depend on the adopted form factors to account for the off-shellness. This t— pion-exchange 
contribution can be safely neglected at threshold. Above threshold, and because the D and D* 
meson masses are quite similar, the exchanged virtual pion will carry a small energy and therefore 
in good approximation q 2 ~ — q 2 . Having in mind the existence of two p— wave vertices in the 
diagram, each of them proportional to the virtual pion momentum, we expect a large cancellation 
of these with the pion propagator (neglecting the pion mass), which will significantly reduce the 
energy dependence of this contribution. 

From the above discussion, it is clear that the predictions of the model of Ref. 53] far from 
threshold are likely subject to large uncertainties due to the contributions mentioned above (pion 
exchanges in the t— channel, u— cut contributions), and possibly other mechanisms. This is the 
reason, why in this work we have only computed properties of the in medium amplitudes close to 
threshold. 



II. FORMALISM: THE DN AND D*N INTERACTION IN NUCLEAR MATTER 

We will calculate the self-energy of the D and D* mesons in nuclear matter from a self-consistent 
calculation in coupled channels that treats the heavy pseudoscalar and vector mesons on equal 
footing, as required by HQS. To incorporate HQS to the meson-baryon interaction we extend 
the WT meson-baryon lagrangian to the SU(8) spin-flavor symmetry group [5^ |. We start from 
the traditional three flavor WT Lagrangian, which is not just SU(3) symmetric but also chiral 
(SUl(3) (g> SUr(3)) invariant. Symbolically, up to an overall constant, the WT interaction is 

C WT = Tr([M\M][B\B}), (1) 



2 For instance, the t— channel p exchange cut for the elastic Z)*E* amplitude goes from about 0.53 GeV to 4.45 GeV. 
Obviously, one does not need to consider explicitly this cut when studying the elastic irN scattering in the region 
of the A(1232) resonance. 

3 Previous models [H, [35, 36] also ignored this contribution because in these works D* -degrees of freedom were not 
explicitly taken into account in the external legs. But, as we have commented, ignoring such degrees of freedom 
was totally unjustified. 
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where mesons (Af ) and baryons (B) fall in the SU(3) representation 8, which is the adjoint represen- 
tation. The commutator indicates a t— channel coupling to the 8 a (antisymmetric) representation. 
For the SU(8) spin-flavor symmetry, the mesons Af fall now in the 63 (adjoint representation) and 
the baryons B are found in the 120, which is fully symmetric. The group reductions 

63 ® 63 = 1 © 63 s © 63 a © 720 © 945 © 945* © 1232 
120 ©120* = 1© 63 © 1232 ©13104 (2) 

lead to a total of four different t— channel SU(8) singlet couplings, that can be used to construct 
s-wave meson-baryon interactions 

((Aft © M)i ® (fit 8) , ((Aft ® Af ) 63a {B ] © 5)63) x , 

((Aft ® M)as. ® (B f ® B) e8 ) , ((Aft © M) 1232 ® (-Bt ® 5)1232) • (3) 

To ensure that the SU(8) amplitudes will reduce to those deduced from the SU(3) WT La- 
grangian in the (8i)meson-(82)baryon subspace (denoting the SU(3) multiplets of dimensionality 
n and spin J by n2j+i), we set all the couplings in Eq. ^ to be zero except for 

£^ 8) = ((Aft ® M) 63o ® (fit ® B) 63 ) 1 , (4) 

which is the natural and unique SU(8) extension of the usual SU(3) WT Lagrangian. To compute 
the matrix elements of the SU(8) WT interaction, £^^ 8 \ we use quark model constructions of 
hadrons with field theoretical methods to express everything in tensor representations as described 
in Appendix A of Ref. [57]]. Thus, we get the tree level amplitudes (we use the convention V = —£)'■ 



where the last factor is due to the spinor normalization convention: uu = vv = 1, as in Refs. 0,0- 
In the above expression IJSC are the meson-baryon isospin, total angular momentum, strangeness 
and charm quantum numbers, Af (E) the common mass (CM energy) of the baryons placed in the 



120 SU(8) representation, and D IJSC a matrix in the coupled channel space (see Ref. 57]). 

However, the SU(8) spin-flavor is strongly broken in nature. The breaking of SU(8) is twofold. 
On one hand, we take into account mass breaking effects by adopting the physical hadron masses 
in the tree level interactions of Eq. ([5]) and in the evaluation of the kinematical thresholds of 
different channels. On the other hand, we consider the difference between the weak non-charmed 
and charmed pseudoscalar and vector meson decay constants. Then, our tree level amplitudes now 
read 



v usc ( r s usc ly/s-Mg-Mt E a + M a E b + M b 

where Af a (Af,) and E a (E b ) are, respectively, the mass and the CM energy of the baryon in the 
a (b) channel. We focus in the non-strange (S = 0) and singly charmed (C = 1) sector, where the 
DN and D*N are embedded. In particular, we look at / = and 1 = 1 channels for J = 1/2 
and J = 3/2. The channels involved in the coupled-channel calculation are given in Table [H where 
below every channel we indicate its mass threshold, Af + m, in MeV units. Compared to [57]], 
we take toa = 1232 MeV, instead of the pole mass. As a consequence, resonances that couple 
strongly to channels with A component might slightly change its position and width, as the case of 
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TABLE I: For each isospin I J sector, all the involved baryon-meson channels are compiled, with their mass 
thresholds, M + m, in MeV shown below. 



1 = 0, J = 1/2 



E c tt ND A c 77 ND* E C K A c w E' C K AD S 

2591.6 2806.15 2833.97 2947.54 2965.11 3069.11 3072.51 3084.18 

AD* E c p A c rj' E> A^4> %~K* WjT* E*K* 

3227.98 3229.05 3244.24 3293.46 3305.92 3361.11 3468.51 3538.01 

I = 1, J = 1/2 



A c tt E c tt ATD ND* E C K E c ?y A c p S^AT T,D S E c p E c w 

2424.5 2591.6 2806.15 2947.54 2965.12 3001.07 3061.95 3072.52 3161.64 3229.05 3236.21 

AD* E> E> YJD* E C K* E^r/ E' C K* E^ E*£>* E^ S*iif* 

3240.62 3293.46 3300.62 3305.45 3361.11 3411.34 3468.51 3473.01 3496.87 3537.42 3538.01 

1 = 0, J = 3/2 



E*tt ND* A c lo S*A AD* E cP E*p A c </> S C A* S' C A* E*#* 
2656.01 2947.54 3069.11 3142.01 3227.98 3229.05 3293.46 3305.92 3361.11 3468.51 3538.01 

I = 1, J = 3/2 



E*tt A\D* A cy o E*r? A£> S*A E c p E c w AD* E*p 

2656.01 2947.54 3061.95 3065.48 3099.23 3142.02 3229.05 3236.21 3240.62 3293.46 

E*w E~DJ E*£> s E C K* E' C K* E^ E*r?' E*£>* E^0 2;! AT* 

3300.62 3305.45 3353.06 3361.11 3468.51 3473.01 3475.75 3496.87 3537.42 3538.01 



(I = 1,J = 1/2) S c (2556) and (J = 1,J = 3/2) S c (2554), few MeV's above the values in Ref. |57(]. 
We also take mx* = 892 MeV. Besides, we also use experimental, when possible, or theoretical 
estimates for the meson decay constants. The used values in this work can be found in Table II of 
Ref. 0. 

With the kernel of the meson-baryon interaction given in Eq. Q, we obtain the coupled DN 
and D*N effective interaction in free space by solving the on-shell Bethe-Salpeter equation 111. Il4l. 

H 0, HI H3 



T "(^= i-y»(^g.W(Vi) y "^- < 7 » 

in the coupled channel space. Here G° is a diagonal matrix consisting of loop functions. 

The free space loop function for channel a reads 

G° a M (VI) = i2M a / Dl (P _ q) D ° Ma {q) , ( 8) 

D° Ba (P-q) = {(P-qf-Ml + iey 1 , (9) 
D°mM = tf-ml+ie)- 1 , (10) 

where s = P 2 , D° is a free hadron propagator, a runs for the allowed baryon-meson channels for 
the given I J sector, and M a and m a are the masses of the baryon B a and meson A4 a in the channel 
a, respectively. The previously defined loop function is ultraviolet (UV) divergent. However, the 
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difference G°(y / si) — G Q {- s /s~2) is finite for any finite values of s± and 82, hence, the function can 
be regularized by setting a finite value of G° at a given point. We choose 

G J IJ \V~s = ^ a J ) = 0, (11) 



with index a running in the coupled channel space, as done in Refs. 36j, 1371 . l57{ |. In those works, 
the subtraction point was taken to be independent of a and J as 



fi 1 ) 2 = a(m 2 th + Ml) , (12) 



where and M t h are the meson and baryon masses of the hadronic channel with lowest mass 
threshold for a fixed I and arbitrary J. The value of a = 0.9698 was adjusted to reproduce 
the position of the well established A c (2595) resonance with IJ = (0,1/2). Then, the same 
value will be used in all other sectors. In this work we will follow the same prescription, taking 

(^=0)2 = a ( MEc 2 + m 2) and (/i /=l)2 = a(MAc 2 + m 2y 

Hence, the renormalized (finite) loop function finally reads: 

G°(V5) = i2M J (D%(P - q) D° M (q) - D B (P - q) D° M (q)) , (13) 

with P and P defined such that P 2 = s, P 2 = (p 1 ) 2 , where, for simplicity, the obvious isospin i", 
spin J and channel a labels have been omitted. This is the (standard) method we use to renormalize 
the free baryon-meson loop. 

The properties of D and D* mesons in nuclear matter are obtained by incorporating the cor- 
responding medium modifications in the effective DN and D*N interactions. One of the sources 
of density dependence comes from the Pauli principle acting on the nucleons. Another source is 
related to the change of the properties of mesons and baryons in the intermediate states due to 
the interaction with nucleons of the Fermi sea. 

Those changes are implemented by using the in-medium hadron propagators instead of the 
corresponding free ones. Therefore, we should define a consistent renormalization scheme, similar 
to that adopted in the free case, for the loop function in a nuclear medium with density p. 

Let us first review what has been done before. The in-medium loop function, G A , used in 
Refs. Isl, 4^|, depends on a regularization cutoff A that renders the UV divergence finite, and it is 



defined as: 



Gft(P) = i2M jT D^{P - q) D P M {q) , (14) 



where Dg, Dj^ are the hadron propagators in a medium with density p. In those works 38l. 1491] . 
the three-momentum cutoff A is fixed in such a way that the free (p = 0) results reproduce certain 
known experimental results (for instance, the position of the resonance A c (2595) for the sector with 
quantum numbers IJSC = 0,1/2,0,1). This way of regularizing the UV loop function induces 
medium corrections of the type: 

G" K {P) = G° A (V~s)+5G p A (P), 

5G A (P) = G p A (P) - G° A (V~s) = i2M jf (D P B {P - q) D P M (q) - D° B (P - q) D%{q)) (.15) 

In this work, we do want to avoid finite-cutoff effects. So we define the in-medium loop function 
as the free one G°, given in Eq. (|13|) and defined as in Refs. 0,0113, l55l . l56l . |57| | without having 
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to introduce any cutoff, plus a term that accounts for the same kind of medium effects as those 
displayed in Eq. (Tl5j) . but taking A large enough (that is A — > oo). Hence, we will use: 

GP{P) = G°{yfs) + 8G P {P) , 
8G P {P) = }™JG P A (P) = i2M J -^L (Lfi B (P - q) D%(q) - D%{P - q) D%(q)) . (16) 

The UV finite 5G P correction contains all the nuclear medium effects affecting the loop, and it is 
independent of the selected subtracting point used to regularizate it. Then, a cutoff is not needed 
for the calculation. The defined loop function at finite density, G p , can be rewritten as 

G P (P) = \2M J -^L (D^{P - q) D» M {q) - D%(P - q) D%{q)) , (17) 

where the integrand is the difference between two terms. The first one corresponds to the baryon 
and meson propagators calculated at density p and total momentum P. The second term depends 
neither on the density nor on P, and it is constructed out of propagators evaluated at p = and 
with fixed total momentum P, such that P 2 = (p 1 ) 2 - This shows that our prescription of Eq. f) 161) 
for G p amounts to assume that the employed subtraction to make the UV divergent function finite 
is independent of the nuclear density. 

For practical numerical purposes we will calculate G p as specified in Eq. (|16p . where the free part 
is analytical, regularized with a subtracting constant and well known, and the medium modification 
part 5G P is numerically evaluated in the same way as done in previous works for G^, but taking 
into account that it has a subtracting part G° providing a 5G P that is UV finite. 

For the DN and D*N channels, we consider Pauli blocking effects on the nucleons together with 
self-energy insertions of the D and D* mesons. The self-energy is obtained self-consistently from 
the in-medium DN and D*N effective interactions, T p d(d*)n > as we will show in the following. 
The corresponding in-medium single-particle propagators are given by: 



2E N (p) \p° - E N (p) +ie p°- E N (p) - ie -p° - E N (p) + ie 
d n(P) + 2m n xP) 




where E^{p) = \jp 2 + , w(<f) = y 1 q 2 + ^f)^*) > ^d(d*)(.Q° ^ Q ) is the D(D*) meson self- 
energy and Sj)m*\ the corresponding meson spectral function. In a very good approximation 
the spectral function for D can be approximated by the free-space one, viz. by a delta function, 
because for that case C = —1, there are no low-lying baryon resonances. Finally, n(p) is the 
Fermi gas nucleon momentum distribution, given by the step function n(p) = H(1cf ~ \p\)i with 
k F = (3vrV2)) 1/3 . 

Using Eq. (|16f) and performing the energy integral over q°, the DN and D*N loop functions 
read 



G p D(D * )N (P) = G° D{Dt)N (V~s) + J 



d 3 q M N 



(2vr) 3 E N {p) 



n{ P ] + (20) 



(P° - E N {p)f - uj(q) 2 + ie 



P° -E N (p) -to{q) +ie Jo P° - E N (p) - w + is 



p=P-q 
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where the first term of the integral, proportional to —n(p), is just what we call Pauli correction 
and accounts for the case where the Pauli blocking on the nucleon is considered and the meson 
in-medium selfenergy is neglected. The second term, proportional to (1 — n(p)), is exactly zero if 
the meson spectral functions <Sd(d*) were taken to be the free one, S^ d »Juj, q) = 8{u—uj{q))/{2u). 
Then, it accounts for the contribution of the in-medium meson modification to the loop function. 
Note that, compared to the Refs. 38, |4^|, we also include the antiparticle contributions in the 



propagators. 

As for DA and D*A channels, we include the self-energy of the D and D* mesons. Then, the 
equivalent of Eq. (I20p for those channels read 



G p D(D * )A (P) = G° D{D , )A (V~s) 



(21) 



+ 



d 3 q M A 



-V(M?)) 



(2vr) 3 E A {p) \P°- E A (p) - u(q) + is 



+ 



dw 



Sd(d*)(^>,q) 



P° -E A {p) -uj + ie 



For the other channels that couple to DN and D*N (see Table H|), we refrain from including 
any medium modifications in the loop function and, therefore, we use the free-space one given in 
Eq. (|13p . This is due to the lack of knowledge on how the properties of some mesons, such as p 
or lo, change in the medium. Only pions in nuclear matter have been intensively studied 6Q. l6l|. 
but, as indicated in Ref. [53] and discussed in the next section, the coupling to intermediate states 
with pions is of minor importance for the dynamical generation of the baryon resonances in the 
5 = and C = 1 sector that governs the DN and D*N dynamics in nuclear matter. 

We can now solve the on-shell Bethe-Salpeter equation in nuclear matter for the in-medium 
amplitudes 



1 



V IJ (V~s) 



l-V IJ (^) GpV j \P) 



(22) 



The in- medium D and D* self-energies are finally obtained by integrating T p ' d(d*)n over the 
nucleon Fermi sea, 



d 3 p 



(1=0,7=1/2)^0 
DN 



(P U ,P) + 3T^ 



(/=l,J=l/2)/ p0 p 



P»,P) 



U D *(q°,q) =f$Ln(p) 



T^ J=1/2) (P°, P) + T^> J=1/2) (P°, P) 



+ 3 D*N 



(1=0, 7=3/2), D o 



P U ,P) + 2T% 



(7=1,7=3/2), d0 



JV 



(23) 



(24) 



where P° = q° + E]y(p) and P = q + p are the total energy and momentum of the DN (D*N) pair 
in the nuclear matter rest frame and the values (g°, q ) stand for the energy and momentum of the 
D and D* meson also in this frame. The Rd(d*)(q°>q) has to be determine self-consistently since 
it is obtained from the in-medium amplitude T^ D ,^ N which contains the D{D*)N loop function 

*-*D(D*)N' an< ^ ^ ms ^ as * Q uan tity itself is a function of ^-d(d*){q°^ q)- From this we obtain the 
corresponding spectral function to complete the integral for the loop function G p d(d*)(n,A)(P° > P ) 
as given in Eqs. (|2QI2ip . 



III. RESULTS 

We start this section by displaying in Fig. Q] the squared amplitude of the D*N — > D*N 
transition for different partial waves as a function of the center-of-mass energy P° for a total 
momentum |P| = 0. In particular, we show certain partial waves and energy regimes where we 
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FIG. 1: Squared of the D*N — > D*N amplitude for different partial waves as function of the center-of-mass 
energy P° for fixed total momentum \P\ = 0. Several resonances are shown: (I = 0, J = 1/2) A c (2595), 
(I = 1,J = 1/2) £ c (2823) and S c (2868), (I = 0, J = 3/2) A c (2660), (I = 0, J = 3/2) A c (2941), (I = 1, 
J = 3/2) £ c (2554) and {1=1, J = 3/2) S c (2902). 



can find seven resonances predicted by the SU(8) model [57j that have or can have experimental 
confirmation [gj, i.e., {I = 0, J = 1/2) A c (2595), {I = 1,J = 1/2) E c (2823) and E c (2868), (I = 0, 
J = 3/2) A c (2660),(7 = 0, J = 3/2) A c (2941), (J = 1, J = 3/2) S c (2554) and (J = 1, J = 3/2) 
S c (2902) resonances. All of them couple to the D*N despite not being the dominant one for 
A c (2660), S c (2823) and S c (2554), as discussed in Ref. [57j]. However, we choose to display these 
amplitudes for different nuclear densities, since they the determine the D* self-energy, as follows 
from Eq. (|24p . We analyze three different cases: (i) solution of the on-shell Bethe-Salpeter equation 
in free space (dotted lines), which was already studied in Ref. [53], (ii) in-medium calculation of 
the on-shell Bethe-Salpeter including Pauli blocking on the nucleon intermediate states at normal 
nuclear matter density po = 0.17 fm -3 (dashed lines), (iii) in-medium solution which incorporates 
Pauli blocking effects and the D and D* self-energies in a self-consistent manner for three densities, 
ranging from 0.5 to 2 po (solid lines). 

The A c (2595) resonance is predominantly a D*N bound state in contrast to SU(4) TVME 



18, m, 37 



48|, |43| . Pauli blocking 



model, where it emerged clS Si DN quasi-bound state 
effects on the intermediate nucleon states move the resonance to higher energies, as already found 
in previous in-medium models 35j, | 



481 ] . due to the restriction of available phase space in the 



unitarization procedure. The shift in mass is of the order of 12 MeV, while for the SU(4) TVME 



model of Ref. 38f] the shift is around a factor of two larger. This change in the energy shift 



provided by Pauli blocking can be attributed either to the different SU(4) and SU(8) kernels, or 
to the renormalization scheme employed to make finite the loop function or the combination of 
both effects. In order to disentangle between them, we study in Fig. [2] the A c (2595) resonance in 
the DN -> DN transition by using the SU(4) TVME model (with S DA r=0) of Ref. [381 for cases 
(i) to (iii) above. We compare results obtained by using the cutoff regularization [38j] and our 
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FIG. 2: Squared of the DN -> DN amplitude in (I = 0, J = 1/2) A c (2595) channel as function of the 
center-of-mass energy P° with \P\ =0. We show results obtained with the SU(4) TVME model with a 
three-momentum cutoff (A = 787MeV/c), as in Ref. [38J, and for comparison, we also show results using 
the renormalization scheme assumed in this work (a = 0.895) [Eqs. (|16I17|) ] . We examined three different 
cases: free amplitude, calculation including Pauli blocking effects at p = po and in-medium self-consistent 
solution at po- 



new renormalization scheme [Eqs. (|16ll7p ] 4 . A similar mass shift for the A c (2595) is observed in 
both SU(4) and SU(8) models when Pauli blocking effects are included for the new renormalization 
scheme. Therefore, we conclude that the different renormalization of the loop function is the main 
source of discrepancy, i.e., the in-medium solution depends strongly on the correct treatment of 
the -y/s dependence of the loop function. The treatment proposed here clearly improves over those 
based on the use of a cutoff. Thus, the resonances which are far offshell from their dominant 
channel will be heavily affected, as in the case of A c (2595). The self-consistent procedure moves 
the resonance closer to the free position because the repulsive effect of the Pauli blocking is tamed 
by the inclusion of the D and D* self-energies. In this case, the difference between SU(8) and 
SU(4) models is not only due to the renormalization scheme but also to the inclusion of the D* 
self-energy, which compensates the attraction felt by the D mesons, as we will see in the next 
figures. 

The (I = 1,J = 1/2) S c (2823) and S c (2868) are shown in the top right panel of Fig. E 
Although no experimental evidence of those resonances is available yet, they lie very close to the 
DN threshold and, therefore, changes in the nuclear medium will have an important effect on 
the D self-energy, which is a matter of interest in this work. In fact, those resonant states are 
significantly modified in the medium, since both resonances couple significantly to DN and D*N 
systems as well as D*A and any medium modification in those systems alters their behavior. 

The next resonance predicted by the SU(8) model [57|, the (/ = 0,J = 3/2) A c (2660), might 
be identified with the experimental A c (2625), which is the charm counterpart of the A(1520). This 



We fix now a = 0.895 to obtain the correct position of A c (2595) 
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state couples strongly to the £*7r channel and more weakly to the D*N pair. Pauli blocking and 
self-consistency have smaller effects than in the case of the A c (2595) resonance. This is because, 
while the A c (2595) resonance varies in the medium due to the changes affecting its two dominant 
channels, D*N and DN, the A c (2660) is only modified via the secondary D*N channel. We 
do not include any medium modifications affecting to its dominant channel, vrS*. Though the 



pion self-energy [60, [6l|] may induce some changes in this resonance, the expected effect of those 



modifications in the D and D* self-energies are minor. On one hand, this partial wave will only have 
a direct contribution to the D* self-energy, and thus the D self-energy will be affected indirectly via 
the simultaneous self-consistent calculation of the D and D* self-energies. On the other hand, the 
effect of this resonance in the D* self-energy is marginal because it only reflects in the low-energy 
tail, far from the quasiparticle peak. Then, we refrain to introduce any medium changes for this 
channel. 

The (/ = 0,J = 3/2) A c (2941) resonance might be a candidate for the experimental A c (2940), 
which J p is unknown [62| . This correspondence is made under the assumption that our model needs 
an additional implementation of p-wave interactions in order to explain the decay into D°p pairs 
reported in Ref. [7J, which is also hinted by the dominant coupling to D*N. This strong coupling 
changes its properties significantly when medium modifications are implemented. Moreover, the 
fact that the A c (2941) lies so close to the D*N threshold will have important consequences on the 
D* self-energy and, hence, on the spectral function, as we will see in the following. 

The (/ = 1,J = 3/2) £ c (2554) resonance has not yet a experimental confirmation. However, 
similarly to the A c (2660), this resonance might be the counterpart in the charm sector of the 
£(1670). It couples strongly to DA and D*A channels, which correspond to the KA in the 
strange sector. Pauli blocking effects on the nucleons are relatively weak because the coupling to 
the D*N channel is approximately half the coupling to the two dominant channels. Changes due 
to the selfconsistent procedure are comparable to the case of the A c (2660). 

The £ c (2902) resonance in the (/ = 1, J = 3/2) sector can be a candidate for the S c (2800) 
resonance, by varying slightly the renormalization scale and if this resonance could be also seen 



in 7T7rA c states 57)]. This is in contrast with SU(4) TVME models, which predict it in the / 



1, J = 1/2 channel 36|, [37J, [38J, |48j]. With regard to medium effects, those are comparable to the 



A c (2595) case. The dominant channel for the generation of this resonance is D*N. Moreover, this 
resonance lies 50 MeV below the D*N threshold. Therefore, modifications due to Pauli blocking 
and self-consistency are expected to be more important than for A c (2660) and £ c (2554), and turn 
out to be comparable to the changes in A c (2595). 

We show next in the left and right panels of Fig. [3] the D and D* self-energies, respectively, 
as functions of the meson energy q°. The D and D* self-energies result from the integration over 
the DN and D*N amplitudes, respectively, after self-consistency is reached simultaneously. In the 
upper panels we display the real part of the self-energies for po at meson zero momentum (solid 
lines) together with the partial wave decomposition (dashed and dash-dotted lines). The partial 
waves weighted by the corresponding factors in Eqs. (I23p and fj24f) are summed up in order to 
obtain the total self-energies. In the lower panels we show the imaginary part of the self-energies 
for densities ranging from 0.5 po to 2 po and q = MeV/c (solid lines). The dotted vertical lines 
in the upper panels indicate the free D and D* meson masses. 

With regard to the D meson self-energy, we observe that in the I = 0, J = 1/2 partial wave 
the contribution of the A c (2595) resonance clearly appears for energies of the D meson around 
1650 MeV, while the resonant state S c (2556) governs the / = 1, J = 1/2 partial wave around 
q° = 1615 MeV. This state couples mostly to D*A, mixed with DN and D*£*, and it is absent 



in the SU(4) models [181 . l36l . l37l . 1381 . 1 481 ]. which do not include channels with a vector meson and 



a 3/2+ baryon. Close to the DN threshold (y/s = 2806 MeV in free space), the / = 1, J = 1/2 
is the dominant partial wave. This is a consequence of the fact that this threshold lies very close 
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FIG. 3: Real and imaginary parts of the D and D* self-energies as functions of the meson energy q°, 
including the decomposition in partial waves (first row), and for different densities at q = OMeV/c (second 
row). The positions of the D and D* meson masses are also shown as vertical lines for reference. 



to a resonant state in the sector I = 1, J = 1/2 of 2823 MeV with a width of V = 35 MeV. 
This resonance is affected by Pauli blocking and self-consistency, i.e, by the nuclear medium, as 
it couples strongly to states with D, D* and nucleon content. Closer to that structure, we also 
found a very narrow state in / = 0, J = 1/2 with a mass of 2821 MeV. This resonant state is 
less modified in the medium, since it couples marginally to the DN channel. Due to its narrow 
width even in nuclear matter, the main contribution to the D meson self-energy close to the DN 
threshold comes from the J = 1/2 S c (2823) resonance but modified by the near resonance J = 1/2 
£ c (2868). Those resonances lie above the DN threshold and, hence, have an attractive effect at 
the DN threshold. 

The imaginary part of the D meson self-energy, in absolute value, grows with increasing density 
because of an enhancement of collision and absorption processes. The change in density can be 
seen more easily in the s pec tral function, as we will show below in Fig. 01 Compared to previous 
results in nuclear matter 35, 38, 47, 48|, 4^], the density dependence of the D meson self-energy are 
qualitatively similar. However, in the SU(8) model, we have a richer spectrum of resonant states 
which is reflected in the self-energy. While the A c (2595)A r_1 and T> C (2800)N~ 1 determine the D 
meson self-energy in SU(4) models [Isl, |48|, those contributions together with few other resonant- 
hole states around q° = 1860-2060 MeV, such as £ C (2823)./V _1 and S c (2868)A r " 1 , clearly manifest 
also in the D self-energy using the SU(8) interaction [57l ]. 

A novelty of the SU(8) model is that it allows to simultaneously obtain the D and D* self- 
energies. The D* self-energy comes from the contribution not only from the J = 1/2 partial 
waves but also from the J = 3/2 ones of the D*N amplitude. As expected in the J = 1/2 
sector, we find the A C (2595)A^~ 1 and £ C (2556)A^~ 1 components for q° = 1650 MeV and q° = 
1615 MeV, respectively. For higher energies around q° = 1880 - 1930 we sec the £ C (2823)N~ 1 and 
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FIG. 4: D and D* spectral functions as function of the meson energy q° for different densities at q — MeV/c. 
The positions of the D and D* meson free masses are also shown for reference (dotted vertical lines) . 



S c (2868)iV- 1 contributions. In the J = 3/2 sector, we find S c (2554)iV- 1 for q° = 1610 MeV. Close 
to the D*N threshold (y^s = 2947 MeV in free space), around q° = 1960 MeV, the £ C (2902)A^ 1 
excitation becomes dominant. This resonant-hole state mixes with A C (2941)7V~ 1 . The combination 
of both J = 3/2 resonances becomes the dominant contribution close to D*N threshold and has 
a repulsive effect in the D* self-energy, as observed in the upper right panel. Density effects have 
a similar outcome as in the case of the D self-energy. Those effects are better visualized with the 
spectral function. 

The spectral functions for D and D* mesons as function of the meson energy q° are displayed in 
the left and right panels of Fig. 01 respectively. The solid lines correspond to the spectral functions 
at zero momentum from 0.5 po to 2 po fo r \o\ = MeV/c. A dotted vertical line indicating the free 
D and D* meson masses is also drawn for reference. 

The quasiparticle peak of the spectral function, which is defined as 



(qf = Q 2 +m 2 + ReU(io qp (q),q) , (25) 



qp 



moves to lower energies with respect to the free mass position for the D meson as density increases. 
As mentioned in Fig. [3l the presence of the J = 1/2 £ c (2823) and S c (2868) resonances above 
threshold have an attractive effect at the DN threshold. In fact, those resonances can be clearly 
seen on the right-hand side of the quasiparticle peak. In the low-energy tail of the D spectral 
function, for energies around 1600 — 1650 MeV, we observe the S c (2556)A r ~ 1 and A,(2595)A^- 1 



excitations. Moreover, other wider resonances are generated in the SU(8) model 57|and they 



combine to give the total D meson spectral function. In the SU(4) TVME models [38|, |48j, [49J , the 
J = 1/2 £ c (2800)Af -1 fully mixes with the quasiparticle peak while the A c (2595)A r_1 appears at 
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TABLE II: DN and D*N scattering lengths (fm) 







D*N 


J= 1/2 1 = 


0.001 + i 0.002 


-0.44 + i 0.19 


(Born approx.) 


( 0.59 + i ) 


( 1.82 + i ) 


J = 1/2 J = 1 


0.33 + i 0.05 


-0.36 + i 0.18 


(Born approx.) 


( 0.20 + i ) 


( 0.07 + i ) 


J = 3/2 1 = 




-1.93 + i 0.19 


(Born approx.) 




( + i ) 


J = 3/2 1= 1 




-0.57 + i 0.15 


(Born approx.) 




( 0.27 + i ) 



the same energies as in our SU(8) model, as expected. 

The quasiparticle peak of the D* spectral function moves to higher energies with density and 
fully mixes with the sub-threshold J = 3/2 A c (2941) resonance. In the left-hand side of the peak 
we observe the mixing of J = 1/2 S c (2868)A r ~ 1 and J = 3/2 S C (2902)A^~ 1 excitations. Other 
dynamically-generated particle-hole states appear for higher and lower energies, such as J = 3/2 
£ C (2554)A^ 1 . 

Density effects result in a broadening of the spectral functions as the collisional and absorption 
processes increase together with a dilution of the resonant-hole states. This outcome is qualitatively 
similar to previous models for the D meson spectral function 35, 3^, 47, 48. \4$\. 



As already mentioned in one of the above references [49], the low-energy tail of the D meson 
spectral function due to resonant-hole states (l^iV -1 ) might help to understand the J/^> suppres- 
sion in an hadronic scenario. However, it is unlikely that this lower tail extends with sufficient 
strength as far as the J/\I / threshold to explain J/^ suppression only via the DD decay. A more 
plausible hadronic contribution for the J/^f suppression is the reduction of its supply from the 
excited charmonia, \ci0-P) or ty' ', which may find in the medium other competitive decay channels 
[49l | . Such a more broad scenario for the J/ty suppression has been pictured recently in thermal 
models [s^]. On the other hand, the spectral function for the D and D* mesons will influence the 
behavior of dynamically-generated hidden and open charm scalar resonances in nuclear matter, as 



already pointed out in Ref. [64] 



Finally we study the properties of D and D* mesons close to the DN and D*N threshold. 
We first present in Table [II] results for the DN and D*N effective interactions in free space. In 
particular, we give the J = 1/2 and J = 3/2 scattering lengths for 1 = and 1 = 1, 



„IJ _ 1 MN T IJ 

a D(D*)N — ~ 4^~^7!~ D(D*)N^D(D*)N > u ° ! 



at D{D*)N threshold and with Mn the nucleon mass. For the DN effective interaction, we find 
that our 1 = scattering length is negligible compared to the 1 = 1 one, in contrast to the SU(4) 
TVME model of Ref. [48| or the meson-exchange model of the Jiilich group [40]. Moreover, our 
positive scattering lengths indicate the attractive behavior of the D meson self-energy close to 
threshold in contrast with the values of those previous references. The discrepancy with previous 
works has its origin in the different resonant-hole composition of the self-energy close to threshold. 
Moreover, as a new development, we also provide the D*N scattering lengths. We note that the 
dominant repulsive contribution comes from the J = 3/2 partial wave. 

Calculated scattering lengths come out radically different from those deduced within the Born 
approximation. This is not surprising because of the strong character of the meson-baryon inter- 
action in the C = 1 sector, and the existence of resonances, which hint non-perturbative physics, 
close to the D^N threshold (in some cases they are placed even below it). 
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FIG. 5: D and D* optical potentials as a function of p/po for \q\ = MeV/c. 



We can now define the D and D* optical potentials in the nuclear medium as 

um = H^afi , (27) 

which, at zero momentum, can be identified as the in-medium shift of the D and D* meson mass. 
In Fig. the real and imaginary parts of the optical potential are shown as function of the density 
for \q\ = MeV/c. The symbols indicate the calculated values for the optical potentials. We refrain 
from showing finite momentum results for the optical potential due to the uncertainties present in 
our calculation once we move away from threshold as result of the WT interaction together with 
ignoring higher-multipolarity interactions. 

The mass shift at p = po stays attractive for the D meson while becomes repulsive for the D* 
meson with increasing density, in correspondence with the behavior of the quasiparticle peaks in 
Fig. JH Similar values for the D meson potential were obtained for in-medium models of Refs. 



3 a . l3a . 1471 . l48l. 1491] . However, as explained in Fig. HI the origin of the attraction can be traced back 
to a different resonant-hole contribution in the SU(4) models compared to the SU(8) case. 

The dilution with density of the spectral functions give rise to the increase of the imaginary 
part of the optical potential for both D and D* mesons, as this corresponds in module to half of 
the width of the spectral function at the quasiparticle peak. However, the D meson width turns 
out to be much more smaller in the SU(8) scheme than in the SU(4) models. Then, we expect to 
find bound states for the D°-nucleus system [4l[ ■ Looking at the strength of the optical potential 
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for density po, we expect those states to be bound at most by 15 MeV and have half-widths lower 
than 4 MeV. Hence, we expect various states to be observable in different nuclei. It can be of 
interest to study the D°-nucleus spectrum, binding energies and widths predicted by the optical 
potential obtained here and compare it with predictions of other models. It is not clear that the 
D + -nucleus hadronic attraction will be able to overcome the Coulomb repulsion to provide similar 
bound states, that will be a subject of future research. Experiments to determine such a -D°-nucleus 
bound states will be welcome. 



IV. CONCLUSIONS 

We have studied the properties of D and D* mesons in symmetric nuclear matter within a 
simultaneous self-consistent coupled-channel unitary approach that implements the features of 
heavy-quark symmetry. The corresponding in-medium solution incorporates Pauli blocking effects, 
and the D and D* meson self-energies in a self-consistent manner. In particular, we have analyzed 
the behavior of dynamically-generated baryonic resonances in the nuclear medium in the C = 1 and 
5 = sector within this SU(8) spin-flavor symmetric model and their influence in the self-energy 
and, hence, the spectral function of the D and D* mesons. We have also obtained the D and D* 
scattering lengths, and computed optical potentials for different density regimes. We have finally 
compared our results with previous SU(4) models 38, 40, 4^, 4s|, paying a special attention to the 



renormalization of the intermediate propagators in the medium beyond the usual cutoff scheme. 

The SU(8) model generates a wider spectrum of resonances with C = 1 and 5 = content 
compared to the previous SU(4) models. While the parameters of both SU(4) and SU(8) models 
are fixed by the (/ = 0,J = 1/2) A c (2595) resonance, the incorporation of vectors mesons in the 
SU(8) scheme generates naturally J = 3/2 resonances^such as A c (2660), A c (2941), £ c (2554) and 
S c (2902), which might be identified experimentally [62]. New resonances are also produced for 
J = 1/2, as S c (2823) and S c (2868), while others are not observed due to the different symmetry 
breaking pattern used in both models. The modifications of the mass and width of these resonances 
in the nuclear medium will strongly depend on the coupling to channels with D, D* and nucleon 
content. Moreover, the resonances close to the DN or D*N thresholds change their properties more 
evidently as compared to those far offshell. The improvement in the regularization/renormalization 
procedure of the intermediate propagators in the nuclear medium beyond the usual cutoff method 
has also an important effect on the in- medium changes of the dynamically- generated resonances, 
in particular, for those lying far offshell from their dominant channel, as the case of the A c (2595). 

The self-energy and, hence, the spectral function of the D and D* mesons show then a rich spec- 
trum of resonant-hole states. The D meson quasiparticle peak mixes strongly with X C (2823)V~ 1 
and S c (2868)A r ~ 1 states while the A c (2595)A r ~ 1 is clearly visible in the low-energy tail. The D* 
spectral function incorporates the J = 3/2 resonances, and the £ c (2902)iV- 1 and A C (2941)V- 1 
fully mix with the quasiparticle peak. As density increases, these l^A^ -1 modes tend to smear 
out and the spectral functions broaden as the collisional and absorption processes increase. This 
broadening in dense matter might have important consequences for the dynamically generation 
of scalar resonances with hidden and open charm content [oil as well as for excited charmonium 
states for the experimental conditions expected in the PANDA and CBM experiments at FAIR [3]. 
This latter experimental scenario, however, requires the incorporation of finite temperature effects. 

The behavior with density of the quasiparticle peaks is better visualized with the optical poten- 
tials. The D meson potential stays attractive while the D* meson one is repulsive with increasing 
density up to twice the normal nuclear matter one. The attractive and repulsive character of 
the DN and D*N interactions close to threshold, respectively, was already observed in free space 
via the scattering lengths. In particular, the optical potentials with density do not follow the low- 
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density approximation, as expected from the complicated resonant-hole structure of the self-energy. 
The imaginary part in both cases increases with density. Compared to in-medium SU(4) TVME 
models, we obtain similar values for the D meson real part of the potential but much smaller 
imaginary parts. This result can have important implications for the observation of -D°-nucleus 
bound states. Work along this line is in progress. 

Future work also includes the study of the influence of the A self-energy in the in-medium D 
and D* self-energies as well as the inclusion of the width of the vector mesons in the meson-baryon 
channels. Moreover, finite temperature effects are mandatory for the analysis and interpretation 
of the data in the future CBM heavy-ion experiment at FAIR. 

To end up this work, we would like to stress that it would be important to count with an 
improved model in the vacuum, and future research along these lines would be highly desirable. 
When considering the findings of this work, one should bear in mind that the deficiencies of the 
model in the free space would definitely affect to results presented here for amplitudes embedded 
in cold nuclear matter. 
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